Quantizing Billiards with Arbitrary Trajectories 
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The structure of the semiclassical trace formula can be used to construct a quasi-classical evolution 
operator whose spectrum has a one-to-one correspondence with the semiclassical quantum spectrum. 
We illustrate this for marginally unstable integrable and non-integrable billiards and demonstrate 
■ its utility by quantizing them using arbitrary non-periodic trajectories. 
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Periodic orbits form the skeleton around which generic classical motion is built, and not surprisingly, many dynam- 
ical quantities of interest can be computed quite accurately in terms of these invariant trajectories |1J]. They also form 
the basis of modern semiclassical quantization schemes which express the trace of the energy dependent quantum 
propagator in terms of lengths, stabilities and focal points of periodic orbits B. Thus resonances in open systems 
and the energy eigenspectra of closed systems can be determined semiclassically even when the dynamics is chaotic. 

A similar relationship also exists for the trace of the classical propagator, L , though this is known explicitly only 
when the system is hyperbolic j|. Here L* governs the evolution of densities 

r 

L*o^(x) = J dyJ(x-f*(y))0(y) (1) 

where x = (q,p) and f* refers to the flow in the full phase space. We denote by A n (t) the eigenvalue corresponding 
to an eigenfunction n (x) || such that L*o0„(x) = A„(i)</> ra (x). It then follows that the autocorrelation of any 
ON analytic function, A(x(t)) can be expressed as Y^n a n^n(t) where the coefficients {a n } are determined by A(x) 
and the eigenfunctions {0„(x)}. The classical spectrum, {A n (t)}, is thus crucial to understanding the evolution of 
£>y correlations and the first step towards determining this is to evaluate the trace of L*. For hyperbolic systems |j] 

O ■ where the summation over p refers to all primitive periodic orbits, T p is the time period and J p the Jacobian matrix 
J> , evaluated on the orbit. Note that the semi-group property, L* 1 o L* 2 = L tl+t2 , for continuous time implies that the 
k> ' eigenvalues A n (t) are of the form e Xnt . The poles in the Laplace transform of Tr L* thus occur at A n much in the 
same way as the trace of the energy-dependent quantum propagator has poles at the eigenenergies, {E n }, of the 
system. There are important differences though between the traces of the quantum and classical evolution operators. 
The quantum case involves the square root of |det(l — J p r | and in addition has Maslov phases so that the classical 
spectrum in general bears no relationship to the quantum spectrum 

Our intention here is to construct quasi-classical propagators similar to Eq. (|l|), but based on the quantum trace 
formula so that {A n } is directly related to {E n }. This can be used to our advantage by evolving classical functions 
and extracting the semiclassical spectrum from the peak positions in its power spectrum. Importantly, the recipe does 
not involve periodic orbits which often have a highly non-trivial topological organisation due to which a systematic 
computation is difficult. To keep matters simple however, we shall carry out this exercise for both integrable and 
non-integrable systems that are marginally unstable. Compared to chaotic systems, these are simpler to handle by 
far but nevertheless serve to illustrate the essential idea. 

Quasi-classical propagators have been studied before by Cvitanovic and Vattay [||J^] though from a different stand- 
point. They consider an extended dynamical space to evolve the quasi-classical wavefunction and ultimately derive a 
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trace formula that is distinct from the usual Gutzwiller trace formula for chaotic systems, and, whose corresponding 
Fredholm determinant is entire. Thus, they do not consider the question of quantization using arbitrary trajectories. 
A resemblance to the method employed here can however be found in the approach of Heller and Tomsovic Q in that 
they work in the time domain and use the fourier transform of the correlation function to determine the spectrum, 
{En}. There are basic differences though, primarily concerning the nature of the evolution. While we use here the 
classical propagator and its quasi-classical adaptation, Heller and Tomsovic use a semiclassical approximation of 
the quantum propagator. As a consequence, the quantities that need to be computed differ. 

We shall, without much loss of generality, confine ourselves to billiards |l(J . The ones that are marginally unstable 
include the integrable circle billiard and rational polygonal billiards. An important characteristic of these systems is 
the existence of two independent constants of motion due to which the invariant surface is two dimensional. In case 
of integrable motion, this is topologically equivalent to a torus while in the general case, the invariant surface has the 
topology of a sphere with g holes. Examples of integrable billiards are few and generic rational polygonal billiards 
are pseudo- integrable with g > 2. 

Before considering the question of quasi-classical quantization, we first introduce the appropriate classical evolution 
operator. For integrable systems, this is easily defined as 




where 9\ and 62 are the angular coordinates on the torus and evolve in time as 9\ = Wi(I\,l2)t + L with = 
dH(Ii, I-2)/dIi and 1 L = <f> r p.dq. Here T;, i = 1, 2 refer to the two irreducible circuits on the torus and p is the 
momentum conjugate to the coordinate q. 

It is easy to see that the eigenfunctions {<j> n {Q\, 82)} are such that (f> n (6i, 6\) — K n {t)(f) n {6x, 62) where A n (t) — e m ™*. 
On demanding that cV^i, #2) be a single valued function of (#1,6*2), it follows that <^ n (0i,02) = e^" 1 61+712 62 •* where 
n = (711,712) is a point on the integer lattice. Thus the eigenvalue, A n (t) = ex.p{it(niuji + tt^u^)}- We shall evaluate 
this explicitly for the rectangular billiard while discussing the trace of L* and we now introduce the appropriate 
evolution operator that applies also to rational polygonal billiards which are generically non-integrable. 

For both integrable and pseudo-integrable polygonal billiards, the dynamics in phase space can be viewed in a singly 
connected region by executing 2g cuts in the invariant surface and identifying edges appropriately. 




FIG. 1. The singly connected region for an L-shaped billiard consists of four copies with edges appropriately identified. A 
trajectory originating near the 37r/2 vertex in 1 is plotted in configuration space using bold lines and the corresponding unfolded 
trajectory is also shown. The latter consists of parallel segments and the trajectory can be parametrized by the angle ip that 
it makes for example with the qi axis. 
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At a given energy, the motion is parametrised by the angle, ip, that a trajectory makes with respect to one of the edges. 
As a trivial example, consider the rectangular billiard. The singly connected region is a larger rectangle consisting of 
four copies corresponding to the four directions that a trajectory can have and these can be glued appropriately to 
form a torus As a non-trivial example (see Fig 1), consider the L-shaped billiard which is pseudo-integrable with 
its invariant surface having, g = 2. Alternately, the surface can be represented by a singly connected region in the 
plane and consists of four copies corresponding to the four possible directions an orbit can have and these are glued 
appropriately. A trajectory in phase space thus consists of parallel segments at an angle <p measured for example with 
respect to one of the sides. It will be useful to note at this point that the same trajectory can also be represented by 
parallel segments at angles tt — p, tt + p and 2tt — p. In general, the number of directions for representing a trajectory 
equals the number of copies, N, that constitute the invariant surface. 

The classical propagator on an invariant surface parametrised by p is thus 

L\p)o^ q ) = J dq' S(q - q n (p)) 0(q') (4) 

where q refers to the position in the singly connected region and q (</?) is the time evolution parametrised by p as 
described above. We denote by {A„(i; p)}, the eigenvalues of L*(p) and from its multiplicative nature, it follows that 
A n (t;p)=e x ^ t . 

The trace of L takes into account all possible invariant surfaces that exist and hence involves an additional 
integration over p jl3) . Thus 

Tr L l — J dp L\ip) = J dpJ2 eXn{,f)t = J dp J *l%-q*(v)) ( 5 ) 

Clearly the only orbits that contribute are the ones that are periodic. Further, the q integrations are simpler if we 
transform to a local coordinate system with one component parallel to the trajectory and the other perpendicular. 
Thus <5|| (<7|| — gjj) = ^S(t — rT p ) where v is the velocity, T p is the period of the orbit and r is the repetition 

number. Similarly, for an orbit of period rT p parametrised by the angle p p , <5j_(<7j_ — Q± P ) = i Pp)/\dq±/dp\ v=lp 
where \dq±/dp\ v=ip — rl p for marginally unstable billiards. Putting these results together and noting that each 
periodic orbit occurs in general at N different values of p, we finally have 
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Tr L l = N UUw^l-rZ,,) (6) 

p r—l ^ 

where I = tv and the summation over p refers to all primitive periodic orbit families with length l p and occupying an 
area a p . We thus have the analogue of Eq. (||) for integrable and polygonal billiards. 

In some cases, it is possible to interpret the periodic orbit sum in Eq. (^|) starting with the appropriate quantum 
trace formula fl5f| 

^ oo 

53 6{E - E n ) = d av (E) + -j= ^ 53 cos(fcrip - j - vrn p - rv^-) (7) 



Here d av (E) refers to the average density of quantal eigenstates, k = \/E, l p is the length of a primitive periodic orbit 
family, n p the number of bounces that it suffers at the boundary and v p the number of caustics encountered by the 
orbit. Note that in the Neumann case, n p = since there is no phase loss on reflection while for polygonal billiards, 
v p = 0. For convenience, we have set fi = 1 and the mass m = 1/2. Starting with the function 

53 f(VE~l)e-P E " = dE f(^El)e'^ E S(E ~ #») (8) 

n ^ e n 

where f(x) = cos(x — 7r/4) and < e < E , it is possible to show using Eq. (Q) that for polygonal billiards |ll| 

OG t ~\ \TJ\ 

Y^ M-J P S{1 _ Hp) = 2^0 + 2^53/(^0 (9) 



p r=1 Vly/f% 

for 0+. In the above, b = J2 P J2 r Jq dE f {VEl) f (VErl p ) and is a constant It follows from 

Eqns. and (||) and the fact that h(x)S(x — a) = h(a)S(x — a) that 
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Tr 1} = 2irNbo + 2irN^2f(y/E^l). (10) 



where {E n } are the Neumann eigenvalues of the system. Thus X n (cp) = i\JE ri sin(</?) |18|. 

For integrable polygonal billiards, Eq. ( |l0|) can in fact be derived directly from the expression for A„(t) given 
earlier. To illustrate this, we consider a rectangular billiard where the hamiltonian expressed in terms of the actions, 
Ii, I2 is H(Ii,I 2 ) = it 2 {I 2 /L\ + Z|/i|) where L\,L 2 are the lengths of the two sides. With I\ = \[EL\ cos(^)/7r and 
I2 = VEL2 sm((p)/ir, it is easy to see that at a given energy, E, each torus is parametrised by a particular value of 
(p. Thus A n (t;(p) = e ^tVE(n 1 cos( v )/L 1 +n 2S in( v )/L 2 ) and the spectrum \g continuous, parametrized by ip. The trace of 
the classical propagator is thus 



Tr L* = I " d( P e u ^ sin(v+ ^ = 2tt^ Jq{^/EJ) (11) 



where / = 2t\E, tan(/j n ) = n\L2i '{n^Lx) and E n — Tt 2 {n\/L\ + n\/ L\). On separating out n = (0, 0) from the rest, 
restricting the summation to the first quadrant of the integer lattice and noting that for a rectangle bo = 1/4, Eq. ( |l0| ) 
follows 

The classical evolution operator thus serves to determine the Neumann spectrum in polygonal billiards. We have 
confirmed this numerically for the rectangular and equilateral triangle billiards and an outline of the method used 
can be found towards the end of this paper. Appropriate modifications however need to be made for the Dirichlct 
spectrum and for systems that have caustics (the circle billiard is an example) and the construction of the evolution 
operator is then guided by the nature of the quantum trace formula 

The quasi- classical evolution operator linking the classical and semiclassical eigenvalues can be defined as 

L* c (<p)o0(q) = J dq' 6(q - q^^e"^-" <*>*0(q') (12) 

where v(t) = z^(q — ' (^3)) and n{t) = n(q -t (ty?)) count respectively the number of caustics and (phase altering) reflections 
encountered by the trajectory q - *^) in time t. As before, the multiplicative nature of L l qc (Lp) implies that its spectrum 
is of the form {e A "^)*} and we shall now show that for the quasi-classical operator defined in Eq. (|l2|), {A„} has a 
one-to-one correspondence with the appropriate quantum spectrum. 
The trace of L qc is given by 



Tr L\ c = J d^J2 eX " (tp)t =Jd^Jd<i S(q - q'fo>)) e 



— in(t)7r — iu(t) '-. 



(13) 



and the contributions to it are due to orbits that are periodic. Since they occur in families, the integrand is independent 
of q as shown earlier and we have 



OO r. 

Tr L< c = AT £ £ ^(l - ri„) / d V 6(<p 

p r=l P 



^ ) e - <n ( t J *(v)) 7r - i!y (q '(*>))£ ^4) 



N ~~f~ S ( l ~ rl v )e- ml ^- irv ^ (15) 
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Starting with the function J2 n 9(VE^l) exp(— 0E n ), it follows from Eq. (0) that for /3 — > + , 

Tr L\ c = J2 A n(t) = 2nNC + 2<KN^g{sfE~l) (16) 

n n 

where {E n } now refers to the desired quantum spectrum, g(x) = ij ^exp(— ix + 7r/4) and C is a constant. It follows 
from 

27T Jo 

that 
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X n ((f) = iy/E n saa.{tp/2) (18) 

Eq. ( O ) forms the central result of this paper. 

In order to numerically demonstrate our results, we first note that 

L\M ° 0(q) = ^(g-*(^))e- in(9 " t(v))7r - <1/(9 "* (v)) 5 = ]T c e ^^/2)i^( g ) (19) 

where {</> n } are the eigenfunctions of L* c and {c„} are the expansion coefficients for the function </>(q). Thus on 
evolving i/>(i; ip) = </>(g* (yi))e _m ^' ? ' ' (v>))^~-^{g' [ (¥>))f anc | averaging this over i/j, we can extract the quantum eigenvalues 
from the power spectrum of (tp(t)} where (. . .) denotes averaging over tp. In practice though, we use the function 
(V'(i))e _ ^ t to compensate for the finite data set with (3 > 0. 

The first application of the quasi-classical evolution operator is provided for a circle billiard of unit radius in Fig. (2). 
The function 4>(q) is chosen to be a Gaussian inside a cell centred at q c and outside [po[ . The power spectrum, S(k) 
is displayed in Fig. (2) where the arrows mark the positions of the Dirichlet EBK eigenvalues pl| . 
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FIG. 2. Power spectrum of (ip(t)) for a circle billiard of unit radius. The arrows mark the position of the Dirichlet EBK 
eigenvalues. 

A total of 1000 arbitrary trajectories with different values of angular momentum have been used and the the power 
spectrum is averaged over 4 cells. We have also generated the Neumann spectrum of the circle billiard (this is distinct 
from the Dirichlet spectrum) by setting n(t) = 0. 

Finally, we provide results for a non-integrable triangular billiard with base angles (tt/5, 37t/10) in Fig. (3). The 
function (ip(t)} has been generated using 300 arbitrary trajectories and averaging over 150 cells while the power 
spectrum has been obtained with /? = 0.006. We find convergence to be rapid with the number of trajectories used. 
The width of individual peaks depends on >/p and the choice of (3 is dictated by several factors which include the 
degree of smoothening attained by averaging over different cells and the length of trajectories. In the present context, 
the value of is also influenced by the height of the shortest peak (at 18.5) and is chosen such that it is distinct from 
the noisy baseline. As in the circle billiard plf ], peak positions in the non-integrable case do not coincide with the 
exact quantum spectrum. The average error is about 8 percent of the mean level spacing in the triangular billiard 
and we find this to be of the same magnitude as in the eigenvalues obtained from Bogomolny's (22) transfer operator 
( see Fig. (3) ). 
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FIG. 3. Power spectrum of (ip(t)) for a non-integrable triangular billiard with base angles (7r/5,37r/10). The upper set of 
arrows mark the positions of the exact Dirichlet eigenvalues while the lower set is obtained from Bogomolny's transfer operator. 

The neglect of higher order corrections thus introduces errors as in all semiclassical schemes and it is not immediately 
obvious how effects such as diffraction can be taken into account in the quantization recipe developed here. We wish 
to point out however that direct applications of periodic orbits have met with little success in determining individual 
quantum levels of pseudointegrable systems |p3| and the results provided here are thus of some significance. 

In addition to the cases presented here, we have also generated the semiclassical spectrum of a circle billiard with a 
single magnetic flux line passing through the centre as well as other non-integrable triangular billiards using arbitrary 
trajectories fiof . Also, it is possible to generate excited states by suitably choosing the function </>(q) such that it has 
nodes. We have thus been able to generate the first 100 levels of an arbitrary triangle |lQ|. 

We provide below a summary of our results : 

• (i) we have introduced the classical evolution operator for marginally unstable systems and shown that for 
polygonal billiards, the classical spectrum on an invariant surface labelled by ip is {e lv/ ~®" Sln (<p) 4 } where {E n } is the 
(quantum) Neumann spectrum. 

• (ii) we have also shown how the classical evolution operator can be appropriately modified to construct a quasi- 
classical adaptation whose spectrum has a one-to-one correspondence with the desired semiclassical quantum spectrum 



• this yields a quantization procedure using arbitrary trajectories and its utility is demonstrated by determining 
individual levels in a class of systems where direct applications of periodic orbits have met with little success so far. 
It is a pleasure to acknowledge several useful discussions with Predrag Cvitanovic and Gregor Tanner. 
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